Isolated hybrid normal/superconducting ring in a magnetic flux: from persistent 

current to Josephson current 
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We investigate the ground state current of an isolated hybrid normal/superconducting ring (NS 
ring), threaded by an Aharonov-Bohm magnetic flux. We calculate the excitation spectrum of the 
ring for any values of the lengths of the normal metal and of the superconductor. We describe the 
nonlinear flux dependence of the energy levels above and below the gap edge. Using a harmonics 
expansion for the current, we isolate the contribution due to these nonlinearities and we show that 
it vanishes for large normal segment length djv. The remaining contribution is very easy to evaluate 
from the linearized low energy spectrum. This decomposition allows us to recover in a controlled 
way the current-flux relationships for SNS junctions and for NS rings. We also study the crossover 
from persistent current to Josephson current in a multichannel NS ring at finite temperature. 
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I. INTRODUCTION 

In spite of the great amount of work devoted to 
this problem, a full understanding of the physics of 
persistent currents in normal mesoscopic rings is still 
lacking. On the other hand, the physics of proxim- 
ity effect in hybrid normal-superconducting mesoscopic 
structures has gained renewed interest recently due to 
progress in nanofabrication techniques. What is the in- 
terplay between these two phenomena? In order to ad- 
dress this question, we study a mesoscopic isolated nor- 
mal/superconducting loop (NS loop) made of a normal 
metal of length and a superconductor of length ds 
as depicted in Fig. [I]. This NS ring is mesoscopic in 
the sense that the normal segment is shorter than the 
coherence length L$. As a consequence of phase coher- 
ence, a non dissipative current flows in the ring when 
a magnetic flux <p is applied. In the absence of super- 
conductor, i.e. for the normal ring, this is the so-called 
persistent current which has the periodicity 4>q = h/e. 
When the superconducting segment is longer than its su- 
perconducting coherence length £ Q , this current is anal- 
ogous to the Josephson current in SNS junctions, with 
periodicity h/2e. 

Biittiker and KlapwijkB(BK) showed that the physical 
mechanism responsible for the crossover between these, 
two cases is the tunneling of Andrccv quasiparticulescl 
through the superconducting segment when its length ds 
becomes of order £ Q . However, BK studied only the low 
energy spectrum where energy levels vary linearly with 
the flux. They constructed the current-flux relationship 
I(4>) in analogy to that of long SNS junctions. In the 
present paper, we investigate the full NS loop spectrum. 
In particular, we describe carefully nonlinear variations of 
the level positions with the flux which appear both below 
and above the gap edge. Then, we address the question 
of a proper calculation of the current 1(0) which includes 
these new spectral features. The standard difficulty is to 
compute the sum of many single level currents —de/d(p 
which arc of the same order of magnitude and alternate 
in sign. In our approach, each harmonic of the total cur- 




FIG. 1: NS ring composed of a normal segment N and a su- 
perconductor segment S, pierced by an Aharonov-Bohm flux 



rent I(4>) is expressed as the sum of a term proportional 
to a single level current at zero energy —de/d(j)(e = 0) 
plus an integral of the "curvature" —d 2 e/d 2 cf> over the 
whole spectrum. The latter term comes specifically from 
the nonlinearities in the spectrum and we show numeri- 
cally that it is significant only for short normal segment, 



typically when d 
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2£ D . For longer normal segment 
2£ , the current /(</>) can be obtained simply from 
the low energy linearized spectrum. As a byproduct, 
one recovers the result for the long SNS junction when 
ds — > oo and for the normal ring when ds — > 0. This 
derivation of I((j>) for long SNS junctions is surprisingly 
simple compared to the original derivation of Bardeen- 
Johnson, Ishii and Svidsinski et a/.Emti. These authors 
found that /(</>) is triangular, but disagreed with the 
value of the critical current. The simplicity of our ap- 
proach enables us to show that Bardeen and Johnson 
found the correct critical current for long SNS junctions. 
Beside the simplicity of the derivation, our decomposi- 
tion gives us the possibility to control the approximation 
since the 'curvature' term is a correction to these well- 
known results. For djv = 0, we can evaluate both terms 
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and by summation, we reconstruct the short junction re- 
sult: the current is 2irA/4> sin(Ax/2) per channel&Q. As 
a new result, we study the crossover from the persistent 
current to the Josephson current at finite temperature 
for single-channel and for multichannel NS loops. 

The paper is organized as follows: in section we 
recall the expression of the thermodynamic potential in 
terms of the excitation spectrum and we derive the exci- 
tation spectrum of a single-channel NS loop for arbitrary 
djq and dg. In section II], we derive our new decomposi- 
tion of the current and evaluate the contribution of the 



nonlinear flux dependent energy levels. In section IV, we 
apply this approach to demonstrate the BK result valid 
in the case of a long normal segment and we consider the 
evolution of the full spectrum when ds varies. Section [v] 
describes the crossover from long SNS junctions to short 
junctions as d/v is decreased. The contributions to the 
ground state current from levels above and below the gap 



A is discussed in section VI Fi nite temperature effects 
are incorporated in section VII to study the transition 



from a diamagnetic to a paramagnetic behavior at small 
flux when ds is lowered or when T is increased. Finally, 
in section VIII shows how to sum over transverse chan- 



nels to obtain the multichannel case from our study of 
the single-channel case. We conclude in section IX. 



II. EXCITATION SPECTRUM AND 
SUPERCURRENT OF THE NS LOOP 

A. Relation between supercurrent and excitation 
spectrum 

We consider a NS loop of perimeter L made of a super- 
conducting segment of length ds and a normal segment 
of length dN ■ The non-dissipative current flowing in this 
system is obtained by differentiation of the thermody- 
namic potential or Gibbs energy Q(T, /i, <f>) with respect 
to the magnetic flux: 



(1) 



For a system with inhomogeneous superconductivity, 
Bardcen et al. and Bccnakker et al. have shown that 
it is possible to express the thermodynamic potential 
in terms of the excitation spcctrum&El This excitation 
spcctruniris found by solving the Bogoliubov-Dc Gcnncs 
equational3: 



H A(x) 
A*(x) -H* 



u(x) 
v(x) 



u(x) 
v(x) 



(2) 



These equations apply when the normal segment is 
shorter than the phase coherence length L^. For such 
a mesoscopic NS ring, excitations are coherent around 
the whole loop and can be described by electron-like and 
hole-like wavefunctions denoted respectively by u{x) and 
v(x). H Q = (—ihd/dx — qA(x)) 2 /2m + V(x) — /i where 



p is the chemical potential and m is the effective mass 
of electrons and holes common for both superconducting 
and normal parts. A(x) is the vector potential due to 
the Aharonov-Bohm flux, V(x) is the disorder potential 
and x is the coordinate along the loop. In this paper, 
we consider the clean system V(x) — 0. Following BK, 
we choose a "square well" model for the superconduct- 
ing gap: A(x) is zero in the normal region and uniform 
A{x) = A in the superconductor. The thermodynamic 
potential can be written in terms of the excitation ener- 
gies and of the superconducting gap: 

Q{T,ix,4>) = -2T^ deln(2cosh^)p exc (e,4>) 

+ J dx\A{x)\ 2 /g + TrH (3) 

where g is the pairing interaction present in the supercon- 
ducting segment. In front of the first integral, the factor 
2 accounts for spin degeneracy and we choose units with 
kB = 1. In this formula, p e ^ c (e, </>) is the density of ex- 
cited states per spin direction. The last two terms in Eq. 
(||) arc independent of the flux. The first term can be 
interpreted as the Gibbs energy for the semiconductor 
model. In the semiconductor model, states with posi- 
tive energies lie at the quasiparticle energies of the initial 
problem and states below the Fermi level lie at the op- 
posite of the former quasiparticles energies, as depicted 
in Fig. |^. By construction, the spectrum of this semi- 
conductor model is fully symmetric with respect to its 
Fermi level and each state can only be singly occupied. 
For this semiconductor, the flux dependent part of the 
thcrmodynamical potential is: 

Sl(T,ft,<i>)=-Tj deln(l + «*P-£)pM) (4) 

where p(e,(f>) = Pexcd^l^) is the semiconductor density 
of states obtained by symmetrisation of the original den- 
sity of excited states as represented in Fig. |[ 



B. Excitation spectrum for arbitrary NS loop 

1. Derivation of the spectrum 

In this section, we calculate the flux-dependent exci- 
tation spectrum e(<f>) of the NS loop. BK have estab- 
lished this excitation spectrum below the gap A^and for 
a fixed number N of electrons per spin direction!!!. Here, 
we calculate the spectrum below and above A, consid- 
ering as well the cases of a fixed number of particles or 
a fixed chemical potential p = kp/2m. The case of a 
fixed number oLelectrons per spin direction corresponds 
to kpL = 7rATlil. The excitation energies are found by 
solving the Bogoliubov-De Gennes equation (^) for the 
NS loop. The eigenvalue equation is obtained by match- 
ing the two-component wavefunctions and their deriva- 
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FIG. 2: Excitation spectrum (top curve) versus semiconduc- 
tor model spectrum (top+bottom curves) of a large purely 
superconducting ring. 



tives at the NS interfaces, see Appendix X A . In the qua- 
siclassical approximation e <C Ep, the resulting equation 
is: 



7 T ( „ 

cos kpL = r e cos - — =F ^p 
\nvp 



e. 



(5) 



where ip = <j> / <f> S& the reduced flux. The minus sign cor- 
responds to excitations around +kp and the plus sign to 
excitations around — kp . O e is a phase shift due to the 
presence of the superconductor that adds to the phase 
shift associated to the motion in the normal segment. 
The function r c defines an energy dependent rcnormal- 
ization of the Fermi wavevector cosfcp-L — ♦ coskpL/r e . 
The functions r e and e have different expressions below 
and above the ga p. In side the gap e < A, they are given 
by, see appendix [X A : 



r e e 



sinh(2wj e — X e ds ) 



sinh 2irj t 



(G) 



(e + 



iVA 2 - e 2 )/A and A e = 
e 2 /hvp. At zero energy, A = A e= o is the inverse of 
the superconducting coherence length £ G = hvp/ A_ which 



is the characteristic length scale for this systemO. The 
modulus of the complex number (^J) is: 

r e =| cosh A e ds + i cot 2rj e sinh A e (is | (7) 

and the phase 9 e satisfies: 

tan 9 C = cot 2r/ e tanh A e c?s (8) 

Above the gap, Eq. becomes, see appendix XA: 

l0e _ smh(i8k e d s + 2S e ) 
reG ~ sinh 25, 1 j 

= cos Sk e ds + icoth 2<5 sin Sk e ds (10) 



so that O e and r e satisfy: 

r t =| cos Sk e ds + i coth2(5 e sm6k e ds 

and: 



tanG £ = coth2(5 £ tan(Sfc e 'i5 



(11) 



(12) 



where e 2<5 « 



Ve 2 - A 2 )/A and 8k t 



Ve 2 — A 2 /hv F . We choose 6 e to have the same integer 
part as 5k t ds- This implies: 

't'Ax\8k t d s \ , T x ( Sk e d s , l" 

, - . + nlnt \ h - 

tanh 25 f \ tt 2 . 



arctan 



(13) 



The spectrum above the gap e > A was not considered in 
the work of BK. Excitation energies are thus solutions of 
the quantization equation (^|) valid whether the energy is 
above or below A. The complexity is hidden in the energy 
dependence of the functions r e and e given respectively 
by Eqs. (00) and by Eqs. (|Jl^Jl3]). These functions 



are plotted in the Figs. (p^,|l5|) of Appendix |X A| . The 
correspondence between the solutions above and below 
the gap is given by the mapping: 



ITJe 



iSk e 



(14) 



The excitation energies are the positive solutions of Eq. 
(||). They are quantized according to e J (n ± ip), the two 
fonctions e J (y) being: 



dN 



1 arccos 

2tt 2tt 



cos kpL 



(15) 



where j = ±1. All the information about the spectrum 
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FIG. 3: e 3 (y) for j = +1 (solid line) and j = —1 (dashed line) 
in a NS loop with ds = 20£ o , djv = 10£ o , containing an even 
number of electrons per spin direction. 

and its flux dependence is contained in the functions 
e 3 (y). We use these functions to calculate the current 
in the following sections. As an example, we have plot- 
ted in Fig. |3| the two functions e 3 (y) for a NS loop with 
dg = 20£ o and rfjy = 10£ o containing an even number of 
particles per spin direction. 
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2. Linear and nonlinear regions of the NS loop spectrum 



The flux dependent spectrum is obtained by folding 
the curves e J (?/) in the interval [—1/2,1/2]. It is shown 
on Fig. ||for the example ds = 20£ o and d,N = 10£ o . At 
zero flux, the first energy levels e J (n, <p — 0) correspond 
to (n, j) = (0, 1), (1,-1), (1, 1), (2, -1), etc. Now, we ex- 
amine the general structure of the excitation spectrum 
for arbitrary values of <In or ds- We distinguish three 
regions: 




FIG. 4: NS loop spectrum e J ± (n,tp) for j = +1 (thin line) 
and j = — 1 (thick line) ds = 20£o, djv = 10£ o for an even 
(left) and for an odd (right) number of electrons N per spin 
direction. For large energies, the two branches j — ±1 tend 
to coincide. The high energy levels exhibit a parity effect and 
are close to those of the normal ring. The Andreev levels are 
insensitive to the parity of N. 

i) The low energy spectrum e < A: in this limit, £ ~ 
e/A.tanhAds and r e ~ cosh Xds- The resulting form of 
the spectrum equation (^) is: 



cos kpL 



(ad* 
- — =F 27iY> 
Tivf 



(16) 



cosh Xds 

and the energy levels vary linearly with the magnetic flux: 



j , v hv F 
e J ±{n,<p) = — 



j ( cos k F L 

n± <p H arccos — r~r~r~ 

2tt VcoshAds 




0,. 5 



FIG. 5: NS loop spectrum e±{n, ip) in the vicinity of A for 
j — +1 (thin line) and j = —1 (thick line) ds — 20£ o , djv = 
10£o and N even (left) or odd (right). 



where d* = d^ + £ tanhA<i,g is the effective size 
of the normal segment. The low energy flux depen- 
dent spectrum is made of linear sections with slopes 
±evF/d*. There are two sets of levels correspond- 
ing to j = +1 and j = —1. Each set is made of 
equally spaced levels with a common average spacing 
hvF/2d*. The energy shift between the two sets is 
hvp/nd*. arccos(cos k F Lj cosh Xds)- For d s ^> £ c , this 
shift is hvpl^d* . Taking into account the spin degen- 
eracy, there are Ad* /7r£ quasiparticle states inside the 
gap. 

ii) The high energy spectrum e 3> A: then C ~ 
5k e ds — eds/fiVF and r t ~ 1. Eq. (||) takes the asymp- 
totic form: 



cos kpL ~ cos 



- — T27r<^ 
Tivf 



and the spectrum is linear in flux: 



,• hv F 
e J ±{n,cp) = — — 



n ± tp + — arccos(cos k F L) 
2n 



(18) 



(19) 



(17) 



The high energy spectrum has the same linear struc- 
ture as the low energy spectrum but with a smaller 
slope ±evF I L which corresponds to excitations extended 
around the whole loop. Indeed, Eq. ( [l9| ) is the linearized 
excitation spectrum of a purely normal ring of perimeter 
L = ds + di- For a given parity, cos kpL = dbl, the two 
sets of levels j = +1 and j = — 1 are in coincidence, as 
depicted in Fig. |]. 

iii) Between the regions i) and ii), the curvature 
d 2 e/d 2 4> is finite and alternates in sign. It is impossible 
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to linearize the spectrum in this region plotted in Fig. 
[5]. One has to be very careful in evaluating the current 
carried by these levels. This will be done in section [II. 



C. From the normal ring to the SNS junction 

In the limit ds/£ D = 0, we have e = and r e = 1 and 
we recover the quantization condition for a normal ring 
of length L: 



( 

cos kpL = cos - — =F 2?t(p 
\ nvF 

The corresponding linearized spectrum is 



„• hv F 
e±{n,(p) = — 



n ± (p -\ arccos(cos kpL) 

2tt 



(20) 



(21) 



characteristic of a purely normal ringtL3. In Eq. (21), 
±tp stands for excitations with momentum around ±kp. 
The index j = 1 corresponds to hole-like excitations with 
I k |< k F and j = — 1 to electron-like excitations | k |> 
kp. 

In the opposite limit, for large ds/£ OI is vanish- 
ingly small in the gap region and e is simply related to 
the Andreev energy dependent phase shift at a NS in- 
terface by 6 C = 7r/2 — arccos(e/A). Consequently, the 
Andreev level spectrum given by Eq. ( |l5| ) becomes: 



-_{n,<p) 



hvF 
d N 



71 ± If ■ 



1 e 
— arccos — 
2tt A 



J 



1 



(22) 



We can express this latter spectrum with only one quan- 
tum number m = 2n + {j — l)/2, the first excited states 
at zero flux corresponding to m = 0, 1, 2, etc: 

hvp 



e±(m,tp) = 



2d 



N 



m ± 2ip ■ 



1 e 

— arccos — 

7T A 



(23) 



We recognize in Eq. ([23|) the spectrum, discovered by Ku- 
lik for bound states in a SNS junction^, the difference be- 
tween the phases xi and X2 of the superconducting order 
parameters of the leads being xi — X2 = 47ry> = 47T0/0 O - 
This spectrum can be understood in terms of quantiza- 
tion along closed orbits in which one electron propagates 
in one direction along the normal segment, is reflected 
as a hole at the superconductor interface with a phase 
shift xis — arccos(e/A), then the hole goes back along 
the normal segment and is finally reflected as an electron 
with an additional phase shift — X2S ~ arccos(e/A). This 
explains why the level spacing at xis = X2S is hvp /2g?at 
corresponding to a box of size 2d^. This scheme applies 
because the electron cannot tunnel through the super- 
conductor in the regime d$ ^Co- 



lli. HARMONIC EXPANSION OF THE 
SUPERCURRENT 

In the previous section, we have shown that the NS 
spectrum can be linearized for low and for high energies 



compared to A. Close to the gap edge, we have identified 
a complicated nonlinear variation of the levels with the 
flux. Instead of the usual decomposition between current 
carcied-by the Andreev levels and by the levels above the 
gapOELZI, we demonstrate that one can extract a contri- 
bution to the current harmonics specifically due to the 
nonlincaritics, namely a term proportional to d 2 e/d 2 (j>. 
Here, the nonlincaritics occur because a quasi-particule 
experiences an energy dependent phase shift when it goes 
through the NS boundary. Our approach is quite general 
and can be applied for other systems where a no nlinear 
behavior occurs. As an example, in appendix X C, we use 
this formalism to calculate the correction to persistent 
current of the normal ring with a quadratic dispersion 
relation. 



A. Derivation of the main result 

We start within the framework of the semiconductor 
model and we use the functions e J (tp) introduced in the 
previous paragraph. In order to simplify the notations, 
we first write the current for one value of j and omit the 
index j for convenience. Hence, we consider the spec- 
trum e(n±ip) and n = .... —1, 0, 1, .. in the semiconductor 
model representation. We can express the Gibbs energy 
([|) in terms of the double integral of the density of states 
N(e, (j)) defined by: 

7V(e»= f de> f de"p(e"» (24) 
in the following manner: 

0(T,/i,^) = j deN(e,4>) (25) 
leading to the current given by Eq. ([l]): 



de 



4Tcosh e/2T V 



(26) 



The quantities p(e, <j>) and N(e, (j)) are even functions of 
the magnetic flux. Omitting the flux independent part, 
we write the Fourier decomposition of iV(e, <p) as: 



N(e,4>) = \^ N m (e) cos2irrrvp 



m— 1 



(27) 



The density of states of the semiconductor model is given 
by: 

n— 00 

p(e,4>)= 5(e-e(n + <r<p)) (28) 

n— — oc,(7=±l 

Using the Poisson summation formula, one gets the 
Fourier harmonics of the density of states : 



Pm(e) 



dy cos 2-Kmy5(e — e(y)) (29) 
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in terms of e(y), which is given by Eq. ( |l5| ) in the NS 
loop problem. After a double integration, one obtains 
the coefficients N m (e) : 

N m {e) = J V{ \y' ^ my ' d f) (30) 
J_ oa 2ttto dy' 

Finally, the current is given by: 

oo 

I((f>) = I m sin 2irmip (31) 

m— 1 

with: 



Im(T) 



27T771 



dc 



-N m (e) (32) 



<Po J-oo 4Tcosli e/2T 
At T = 0, the current harmonics are given by I m (T 



0) = —2-KmN m (t = 0)/(f> o - Integrating by parts Eq. (30), 
one gets : 



I m (T = 0) 



2 1 



dc 



■Km (j) [ dy 

y ° , d 2 e „ 
dy—^— cos ZTrmy 

-oo dy 



(j/o) cos 2irmy 



(33) 



with y = j/(e = 0). We have assumed that the slope 
de/dy is vanishing at the bottom y = — oo of the semicon- 
ductor valence band. Eq. (33) is the general expression 
of the current for a spectrum e(n ± ip). For the NS loop 
case, we have to sum contributions from the two branches 
of levels j = ±1. To understand Eq. ((33[), we recall that 
y plays the role of the reduced flux ip = <f>/ <p . The first 
term is related to the slope of e(y) at zero energy, i.e. to 
the current —de/dcf> carried by the zero energy Andreev 
level. It leads to a triangular I((f>) current-flux relation- 
ship. The second term is a sum over the whole spectrum 
of an integrand proportional to de/dy, i.e. d 2 e/d 2 (f>. Only 
the region around the gap edge with nonlinearities gives 
a non zero contribution, regardless whether these nonlin- 
earities are located below or above A. For this reason, 
our representation of the current is different from the 
usual decomposition of the Josephson current for SNS 
junctions as a contribution from the discrete spectrum 
below the gap plus a cojitribution from the continuum 
spectrum above the gapEJ. 



B. Numerical evaluation. 




FIG. 6: First harmonic I m =i (top curve) and second har- 
monic 7 m= 2 (bottom curve) in units of A/^„ as a function of 
djv, for different values of ds- Dashed lines represent the BJ 
approximation. Inset shows the difference between the second 
harmonic and its BJ evaluation for ds — 1 and ds = 10. The 
dot represents the expected value for djv =0 and ds ~S> £o- 



first term in Eq. ( p3|) which corresponds to the Bardeen- 
Johnson(BJ) approximation used by BK. One sees that 
the second term is already negligible for djv 2. For 
small values of d^, both terms in Eq. ( |33| ) are of the same 
order of magnitude. The difference between the exact 
result and the BJ approximation, which is shown in the 
inset of Fig. || for I m —2, decreases monotonously with 
increasing d^. We note that this correction to the BJ 
approximation is roughly independent of ds for m = 2. 
As it is expected for odd harmonics, l m —\ decreases as 
ds increases, see Fig. ||. For ds = 10£ o , our numerical 
evaluation of the second harmonic is in good agreement 
with the analytical value 7 TO= 2 = — 16A/30 O expected for 
dpj — and ds = oo from Eq. (fry). 



In this section, we show that the relative weight of the 
two terms in Eq. (33) is related to the value of djy- We 
evaluate numerically the integral term in Eq. ( |33| ) for 
loops with finite lengths < rfjv < 10£ o and < ds < 
10£ o . In Fig. (H), we have plotted the first I m =i and 
second I m —% harmonics of I((j>) as a function of g?at, for 
different values of ds- The dashed lines represent the 



C. Conclusion 

We have identified a term specifically due to nonlineri- 
ties in the spectrum and we have shown numerically that 
it is small if d^ ^ 2£ Q . In the following sections, we use 
Eq. ( |33| ) to recover analytical expressions in the extreme 
cases: d^ 3> £ for any ds (section |Tv| ) and ds £ Q 
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for dpj — > (section |v|). The decomposition (^) is valid 
for arbitrary djq with e 3 (y) given by Eq. (|lq) and can 
be used to calculate the current-flux relationship in any 
NS loop. At finite temperature, we find a crossover from 
paramagnetic to diamagnetic behavior at sma ll flux when 



ds is increased or T lowered, see section VII 



IV. LONG NORMAL SEGMENT 



of a normal ring of perimeter d* and with cos kp L re- 
placed by cos kphj cosh Ads- The low energy spectra for 
ds = 5^ o ,20^ o ,oo are quite similar and correspond to 
quasiparticle motion confined in the normal region of the 
loop. The spectrum above the gap is more sensitive to 
the value of ds because the level spacing scales as 1/L 
corresponding to a quasiparticle motion confined in the 
whole loop of length L = ds + djv- 





. 5 


f = 




0. 


f) 
0.5 


1 


Hv) 




-1 


\l 0.5 



1 /M 



NS loop spectrum (n, <p) for a 



FIG. 7: 

c) ds = 5£ and d) ds = 20£ o , keeping 



d s = 0, b) d s ■■ 
n equal to 10f; o . 



= 6», 

The 

number of electrons per spin direction N is even. Levels with 
j = +1 correspond to the thin lines and those with j = — 1 to 
the thick lines. The currents e), f), g), h) are plotted in units 
of I = 2evF/d* below the corresponding spectrum a), b), c), 
d). 

In this section, we study the evolution of the spectrum 
and of /(</>) as a function of ds for a NS ring with d 



N 



> 



2£ D . We have shown in paragraph III B| that in this case 
the second term in Eq. (p3f) is negligible. 



A. Excitation spectrum for arbitrary ds 

In Fig. 0(a,b,c,d), we have plotted the excitation spec- 
trum for ds = 0, £ D , 5£ , 20£ o keeping djv equal to 10£ o - 
For the low energy part of these spectra, the flux depen- 
dence is linear with the slope ±hvp/dN- For the normal 
ring ds — 0, the branches j = ±1 coincide and we recover 
the spectrum of the normal ring with the double degen- 
eracy due to the spin . For finite ds, these branches 
are shifted by A<p = ± aiccos(cos(kp L)/ cosh. Xds)/ 2n. 
For ds = 5£ and ds — 20£ o , they are shifted by 
half a quantum of flux <j) /2 . Indeed, the low en- 
ergy excitation spectrum is the same as the spectrum 



B. Current-flux relationship for arbitrary ds 

The Fourier coefficients in the harmonics expansion of 
the current Eq. ( |3l| ) are J m = Z+ + 1~ with: 



V 



2 1 d^ 



(y J o) cos 2-nmyl 



(34) 



7rm 4> dy 

and they depend only on the zero energy Andreev level 
The equation (|l5|) yields: 



2-Kyl 



-j arccos 



/ cos UfL 
\ cosh Xds 



(35) 



and de? /dy = hvp/d* = Ea- The energy Ea is the typ- 
ical displacement of one energy level when the flux is 
varied. This is also the energy spacing between Andreev 
levels. The order of magnitude of the critical current is 
then Ea/4> - According to Eq. (j33|), the Fourier expan- 
sion of the current for arbitrary ds reads: 



4 ev F ^ T m (X) . 



7r d* ' — ' m 

m—l 



with 



and: 



X 



cos kpL 
cosh Xds 



(36) 



(37) 



t (X) = cos (m arccos X) (38) 

The coefficients T m (X) are the m-th order Tchcbytchcv 
polynomials. The parameter X depends both on the 
band filling and on the crossover parameter Ads. The 
first Tchebychev polynomials are T\{X) = X, T^X) = 
2X 2 - 1, T 3 (X) = 4A 3 - 3X, .... For a fixed number of 
electrons per spin direction, namely for cos kpL = ±1, 
the beginning of the_current expansion is given by the 
following expressionlij: 



4 evp 
n~d*~ 



±- 



cosh Ads 



sm2Trm(j)/4> 



1 — sinh Ads 



sin 4irm(j) / 4> 



(39) 



2 cosh Ads 

The formula above describes the suppression of the first 
harmonic m = 1 when ds goes to infinity. The suppres- 
sion of odd harmonics is a general feature of the crossover 
from persistent current in normal loops to Josephson cur- 
rent in SNS junctions. 
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C. Case ds 3> £ : Josephson limit 



LONG SUPERCONDUCTING SEGMENT 



We recall that for large ds, the NS loop threaded by a 
magnetic flux <fi behaves like a SNS junction with a su- 
perconducting phase shift 4n<j>/<j> between the leads. In 
this ds — * oo limit, the functions e J (y) are flat outside 
the gap. Due to the infinite size of the system, the spec- 
trum becomes a true continuum above the gap. Far from 
the gap, the density of levels is given by: 



dy 

de 



dN 
hvp 



ds 



hv F (e 2 - A 2 ) 1 / 2 



(40) 



It is obtained from expressions (|l l|) , ( pL3| ) and The 
first term is the density of levels in a normal ring of 
perimeter djv at the Fermi level and the second term 
is the BCS singularity at A. At e ^ A, the total den- 
sity of levels tends to those of a normal loop of perimeter 
L = ds + dj^. Inside the gap for e < A, we have: 



2V» = 



Tivp 



+ (1-3)-. 



(41) 



which corresponds to the Kulik spectrum (|22|). It is dom- 
inated by the linear behaviour of the first term for long 
junctions dN 3> £ , except very close to A. Below the 
gap, the flux dependent spectrum is similar to those plot- 
ted in Figs. f?](c,d). In this limit, X — ► and only even 
m = 2p harmonics are non zero in Eq. ( |36| ) because 
?2p+i(0) = 0. It leads to the following <j) /2 periodic 
current: 



7T dN + t,o 



E 



(-ir 



sm4irp(j)/(f) (42) 



This is the Fourier expansion of a saw-tooth current-flux 
relationship for the purely long itudinal channel of long 
SNS junctions. In section VIII , wc show that it corre- 
sponds to the result of Bardccn and Johnson 0. 



D. Case ds = 0: persistent currents in normal rings 



We consider now the case of large ds and we study 
the crossover from long to short SNS junctions obtained 
by decreasing dN- For vanishing dj\r, the levels acquire 
a nonlinear behavior over a large energy range of order 
A. Then, the two terms in the decomposition Eq. ( |33| ) 
become of the same order of magnitude. For d/v = 0, 
we recover analytically the current through a short SNS 
junction. 



A. Excitation spectrum for arbitrary djv 
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If we remove the superconductor, the parameter X is 
equal to cosUfL and the effective length d* is the perime- 
ter L. Then, the model describes the persistent current 
in a purely normal ring of length L with ffxed chemical 
potential. We recover the well-known resuhxJ: 



FIG. 8: NS loop spectrum e±(n, ip) for ds = 20£ o and a) 
djv = 0, b) djv = £ and c) djv = 10£ o - The currents d) 
and e) are expressed in units of I = 2evF /d* and correspond 
respectively to the spectra a) and c). Levels with j = +1 
correspond to the thin lines and those with j — —1 to the 
thick lines. 



I 00 IT 

4 evp tt-^ cosmkpL 

y Sin 27T771 



7T L 



m— 1 



(43) 



the current 1(4') includes both spin directions. This re- 
sult is correct in the framework of the quasi-classical ap- 
proximation e <C Ep. In fact, Eq. ( |i^ ) is the zero order 
contribution in the 1 /kpL expansion. The following term 
in this expansion is due to the quadratic dispersion rela- 



tion of free electrons and is evaluated in appendix X C 



On Fig. ||(a,b,c), we have plotted the excitation spec- 
trum for dN = 0,£ o , 10£ o and for ds = 20£ o . For dN = 0, 
there is only one spin degenerate Andreev level in the 
gap. As dN is increased, new Andreev levels appear, 
see Fig. B(a) t o Fig. ||(b) and Fig. |8|(c). As discussed 
in section|IBj, the number of Andreev levels is roughly 
4(djv/^o + I)/ 71 "- The spectrum is linear at low energy 
and only the last Andreev level at the vicinity of A varies 
nonlincarly with the magnetic flux. 
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B. Current-flux relationship for djv = 
junctions 



short 



In the short junction limit djy = 0, the eigenvalue equa- 
tion is: 

2iry j (e) = - arccos ± + (44) 

The function e J (j/) is plotted in Fig. There is only one 
Andreev level in the gap which corresponds alternalively 
to j = 1 or j — — 1. The spectrum can be written as: 



e(<f>) = A I cos27T 



(45) 



At T = 0, the current corresponding to this unique An- 
dreev level is 4> /2 periodic and given by: 



A 

-2tt— sin 27r 



(46) 



for \4>/4> \ < 1/4, see Fig. BLd). This result was ob- 
tained by Kulik-OmeryanchukQ and Beenakker and H. 
van HoutenB. In appendix XB, we show how to recover 
this result from our formalism and Eq. (|33|). In this case, 
both terms in this equation contribute with the same or- 
der of magnitude, so that the BJ approximation clearly 
breaks down for g?at = 0. 

Finally, we can compare the short djy case, A < Ea = 
Iivf/cI*, and the long <1n case, A > Ea, developed in 
IV. In both cases, the current is (/) /2-periodic, diamag- 
netic at small flux and there is a jump at </> D /4. In the 
former case, I(4>) is triangular and the critical current 
is of order Ea/4> - In the latter case, I(<p) has the sine 
dependence and the critical current is of order A/<f> . It 
is a particular case of the well-known statement that the 
critical Joscphson current in a SNS junction is giwn by 
the minimum of the two energy scales Ea and AHj, Ea 
being the Thouless energy of the clean NS loop. 




FIG. 9: Contributions to the second harmonic from states 
outside the gap an d from states below the gap I™ = 2 in 

units of A/<t> for ds = £0 and ds = 10£ o as a function of dpj. 
The dot represents the expected cancellation of the current 
carried by states above the gap in the limit of large ds- 



For the case of the short SNS junction, d^ = and 
ds ^ £0, the current is solely carried by the single An- 
dreev level. As seen on Fig. ^, when ds is close to £ Q , 
there is a contribution from states above the gap, which 
is significant when ds is close to £ D . 



VI. WHICH LEVELS CARRY THE CURRENT? 

The equilibrium current in a NS loop is carried boll 
levels below and above the superconducting gap AEBL 
Although our formalism is based on a decomposition be- 
tween linear vs. nonlinear flux dependence of the levels, 
we can separate in each term of Eq. ([33]) into contribu- 
tions from states above and below the gap. The result is 
shown on Fig. ^|. Although the total current decreases 
monotonously with <i/v, both contributions are oscillat- 
ing functions of ci/v. These oscillations in the subgap 
current harmonics correspond to the apparition of new 
Andreev levels below the gap when d^ increases. Since 
the number of Andreev levels is Ad* /tt^q, the oscillations 
have periodicity 7r£o/2. Although we have not checked 
this numerically, we believe that the contribution carried 
by the states outside the gap cancels whenever the level 
crossing the gap has zero slope. Otherwise, when an An- 
dreev level crosses the. gap, it still carries current on the 
other side of the gaplijO. 



VII. TEMPERATURE EFFECT 

We treat the effect of a finite temperature on the 
crossover from SNS junctions to normal rings. We es- 
sentially focus on the case dN ^ 2£o- 



A. Harmonics expansion at finite temperature 



We use Eq. ( p2|) to compute the harmonics of the 
current at finite temperature. Compared to Eq. (|3^), 
these harmonics are reduced by the following thermal 
factor: 



I m {T) = I m (T = 0) 



sinh x 



(47) 



with x = 2ir 2 mT / Ea ■ The characteristic energy scale 
associated with the m-th harmonic is given by Ea/ttl. 
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The resulting current for cIn ^ 2£ Q is: 



C. Ensemble average 



T ^ T m (X) 



, sinli7rm-?^ 

m— 1 A 4o 



■^r sin 27rm0 /0 O (48) 



where X and T m (X) are given by Eqs. (37,3a). 



B. Transition from dia- to paramagnetic loops 

For N even, a normal ring has a paramagnetic mag- 
netization at zero flux while a SNS junction carries a 
diamagnctic current. At T = 0, the transition from dia- 
magnetic to paramagnetic behavior occurs at ds = 0. 
At finite temperature, this transition occurs at a finite 
ds- At fixed ds and c/jv, there is a similar crossover as 
a function of the temperature. In Fig. |l^, we see that 
the slope at the origin of the I(cj>) curve changes sign at 
T = T*(ds,dN). The small flux current is diamagnctic 
for T < T* and paramagnetic for T > T*. The function 
T*(ds,dN) is an increasing function of ds and a decreas- 
ing function of diq given by: 



ird* 



(49) 



This behavior can be understood from the harmonics ex- 
pansion Eq. (f48[). For N even, the first harmonic is 
paramagnetic while the second is diamagnctic. When 
the temperature is increased, the reduction of the second 
harmonic is stronger than for the first one and the result- 
ing current becomes paramagnetic. The I(4>) curve for N 
odd is obtained by a </> D /2 translation of the /(</>) shown 
on Fig. [it]. One sees that the magnetization at small 
flux is always diamagnctic. Indeed, all the harmonics are 
negative for N odd. 



0.25 1(0) 



= 0.04A 



ds = 2?„ 




0.5 

V = <P/<Po 



FIG. 10: Current-flux relationships for ds = 2£ , djv = 10£ o 
and for cos k F L =1. At T = T* = 0.056A, the slope dl/d<f> 
at the origin = changes in sign. The current is represented 
in units of I — A/(f> - 



Here, we calculate the average current of a large as- 
sembly of isolated single-channel NS loops with djv 3> £o 
and a fixed number of particles. Statistically, half of them 
have an even number of electrons per spin direction N 
and for the others N is odd. As we have seen in the 
previous paragraph, a transition from diamagnetism to 
paramagnetism occurs for the rings with even TV, while 
the rings with N odd are always diamagnctic. Does the 
total orbital magnetism of the assembly exhibits a tran- 
sition? We know that an assembly of normal rings has 
paramagnetic magnetization at small flux and is o /2- 
pcriodicEl In the opposite limit ds/£, D 3> lj the current 
in each loop is diamagnetic and 4> /2 periodic. The tran- 
sition occurs at a very small value of ds close to 0.5£ o . 
In conclusion, we find that the current is always diamag- 
netic for values of ds and d_/v above £ Q . 



VIII. MULTICHANNEL RINGS 

Up to now, we have considered a single-channel NS 
ring. From now on, we extend our study to multichannel 
NS rings. First, we present spectra with a small num- 
ber of transverse channels and follow the evolution of the 
different channels when djy and ds are varied. Secondly, 
we study the crossover from Josephson current to persis- 
tent current in the clean multichannel NS loop at finite 
temperature. 

The spectrum of clean multichannel rings can be ob- 
tained straightforwardly since the different channels are 
decoupled and characterized by their momenta k y ,k z 
quantized along transverse directions y and z. The spec- 
trum of each of these channels is simply obtained by 
the substitutions vf — > vp x and kp — > kp x in the Eqs. 

■ As an example, we have represented on Fig. |ll| 
the spectra of a ring with a square section of size (2\fY 
for different values of d^ and ds- As djv — > 0, the 
spectra of the different channels all shrink towards the 
single-channel spectrum and become completely degen- 
erate when dps = 0. As a consequence, the I(<p) for a M- 
channel short junctions djv = is exactly M times the 
single-channel result (46). In long SNS junctions, the 
current I ((f)) carried by any transverse channel has the 
same flux dependence as the longitudinal one, namely a 
triangular I{<j)) with current jumps at <f> = ±0 o /4. Nev- 
ertheless, the critical current is different for each channel 
as we can see from the spectrum Fig. |ll|, because the 
slopes —de/d(j)(e = 0) are different. Finally, the total 
current stays proportional to the number of transverse 
channels M. 

The total current is the sum of the single-channel cur- 
rents. For d N > 2£ D , it is given by: 



E 



ev Fx 



n k~k dN 



to 



— ' m 
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FIG. 11: Multichannel spectra. The cross section is a square 
of size (2Af) 2 - For clarity, we have represented only a few 
flux-dependent channels among the M = k%S/<in transverse 
channels, a) ds = £o and d-N = 10£ o , b) ds = 10£ o and 
djv = 10£ o , c) ds — 10£ o and djv = £ , d) ds = 10£ o and 
d N = 0. 



current from Mevp / '(d>N + Co) for large ds to zero as 
ds approaches zero. This can be understood from the 
spectra Fig. [ll]. Indeed, when ds is large, the low energy 
spectra of each channel are in phase with e = for <p = 
±(/> /4. When ds — ► 0, there is finite dephasing between 
the different spectra given by Aip = l/2n. arccosX which 
leads to a cancellation of the total current. 



h/M 



o . 1 



dsKo 



where: 



X = 



cos kp x L 
cosh Ads 



(51) 



In the case of a NS ring with many transverse channels, 
the discrete sum Eq. ( |50| ) over transverse channels can be 
replaced by an integral. Each channel carries a current 
given by Eq. ( fl8| ) with vp — > vp x = vpcos9, where 
9 is the angle between the direction x and the Fermi 
momentum of the channel. By counting of the (k y ,k z ) 
satisfying k 2 . + k 2 , + k 2 = kp for a given incidence 6, we 
obtain the total current: 



I{<t>) = 16ttM— / d9 cos 2 6 V - 



T m {X) . 

. T j, sin Zirmtp 

sinh ixm-^ — 

(52) 



where M = k F S/A-ir is the number of transverse channels. 
In the limit ds > Co and at T = 0, Eq. (||) leads to: 



AM ev F 
3tt dN + £ ( 



E 



(-i)p 

sin Airpip (53) 

P 



FIG. 12: Critical current per transverse channel as a function 
of d s for d N = 10 and k F £o > 1 at T = 0.03A. 



IX. CONCLUSION 

We have calculated the full excitation spectrum of a 
NS loop threaded by an Aharonov-Bohm flux, for any 
value of div and ds'. in particular, we show for the first 
time the spectrum of the NS loop above the gap. We 
have identified the contribution to the current directly 
originating from the nonlincaritics in the flux dependent 
spectrum. We recover known results for short and long 
SNS junctions when ds Co- For the single-channel 
NS ring at zero temperature, we recover the result of 
BK and our method allows us to clarify and justify the 
approximation made by BK. Moreover, we extend the 
study of the NS loop to the finite temperature and to the 
multichannel cases. 

We would like to thank M. Buttikcr, H. Bouchiat and 
S. Gucron for useful and stimulating discussions. 



The corresponding critical current is proportional to M: 

2M ev F 



3 d 



N 



Co 



(54) 



This_is the old result foruad by Bardeen and JohnsonS. 
Ishiicl and Svidsinski et alu found different numerical pre- 
factors. 

In the opposite limit of the multichannel pormal ring 
ds = 0, the total current averages to zcroE2l for large 
kp^ a . In Fig. 115, we show the the crossover of the critical 



X. APPENDIXES 
A. NS loop excitation spectrum 

We consider a purely one dimensionnal NS loop with 
a superconducting segment in the region < x < ds and 
a normal segment in ds < x < L. This loop is threaded 
by a magnetic Aharonov-Bohm flux </>. 

We first investigate states with energies below A. In 
the normal metal, the quasiparticles wavefunctions are 
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plane waves oscillating at wavevectors k e /h = kp±e/hvF 
just below (hole solution) or just above (electron solu- 
tion) the Fermi momentum +kp: 



D 







ikhX 



In the superconductor , they are: 



u(x) 
v(x) 



C 



(55) 



3 (ifcF+A e )a: 



(56) 



where e 2tr >- = (e + iVA 2 -e 2 )/A and A e 
VA 2 — e 2 /hvF. We express the continuity of these func- 
tions at the NS interfaces. In presence of a reduced flux 
tp = 4>/4>o it leads to the system: 



Ae tk ' L = e 2m *{e"^C + e'" h D) 



fc.de 



ikhds 



,(ikp-\ e )d s e iVc(J _|_ e {ikF+K)<ls e -iv* 
Jjk F -X)d s g-M/e Q _|_ e {ik F +X)d s & ir} t jj 



The continuity of the derivatives is automatically satis- 
fied in the quasi-classical approximation e <C Ep. In this 
approximation, there is no mixing between excitations 
around +kp and excitations around — hp. For this rea- 
son, it was possible to consider only solution oscillating 
around +kp in the ansatz Eqs. ( |55| , |56| ). The determinant 
of the system Eq. (^7|) gives the eigenvalues equation: 

2i sin 2ij e cos kpL = e 7 ™ v_ >kF sinh(Ads + 2ir] e ) ~ c.c. 

(58) 

which is identical to the spectrum obtained by BK: 



cos kpL sin 2r\ € 



sin 2ry e cosh \ e ds cos [ N ~ 2ittp 
\nvF 

— cos 2r] e sinh \ e d$ sin [ \ — — 2ttw 
\nvp 



This formula can be reduced in the following form: 



cos kpL = r t cos [ |— ^- — 2ntp + 6, 
V nvp 



with the parameters r £ and £ defined by: 

e iO E _ sinh(2i?7 e - A e d s ) 
sinh 2i7y € 
= cosh A e ds + i cot 2?y e sinh A e ds 

One expression for the modulus r e is: 

r e =| coshA e (is + i cot 2t7 £ sinh X e ds \ 

and the phase O e is satisfies: 

tan £ = cot 2rj e tanh A e ds 



(59) 

(60) 
(61) 

(62) 

(63) 



We now look for quasiparticle states with energies 
above A. In the normal region, the form of the wavefunc- 
tions is unchanged. In the superconductor, they become: 



v(x) 



i{k F —&k E )x _ 



D 



a i(k F +Sk e )x 



(64) 



where e 25c = (e + s/ e 2 — A 2 )/A. Therefore, the eigen- 
value equation can be obtained directly from the preced- 
ing study with the replacement — ir\ t — > S e and A — > i5k e . 
We obtain the Eq. (59|) with the complex parameter: 




FIG. 13: l/r e as a function of energy for different values of 
ds in units of £ Q . 




FIG. 14: Phase shift O e as a function of energy for different 
values of ds. The curves correspond to ds = 1,2,4,6,8,10 
in units of £ from bottom to top. For ds — > oo, O e — » 
7r/2 — arccos(e/A). 



smh(idk c ds + 25 e ) 



sinh 2S e 

= cos 6k t ds + i coth 26 smSk e ds 

The modulus is given by: 

r c =| cos 5k t ds + i coth 25 sin Sk e ds | 

and the phase by: 

tan0 c = coth 2<5 e tan Sk e ds 

Equations (|67| , |68| ) can be obtained directly from 
by prolongation to energies e > A. The fonctions r e and 
6 C are plotted in Figs. ( |il^ . |i4| , |l5| ) for different values of 
ds- We note that r&e lB& = 1 + idg/^ . We have seen 
in the paragraph [I C that these functions have simple 
limits for very large and for very small ds/£ - 



(65) 
(66) 

(67) 

(68) 
fit 
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FIG. 15: Asymptotic behaviour of O e above A for different 
values of ds- The curves correspond to ds = 1, 2, 4, 6, 8, 10 in 
units of £ from bottom to top. At large energies, e is close 
to 5k e ds which is represented in dashed lines for each value 
of ds . 



Similar calculation for j = — 1, give the same result for 
I m when to > 2. The case m = 1 has to be considered 
separately and it is easy to see that the j = +1 and 
j = —1 cancel each other. Consequently, the current is 
given by: 



4A 



p=l 



■ sin 4irp<fr / (j> (71) 



p 2 - 1/4 

which is nothing but the harmonics expansion of Eq. 
©■ 



Exact calculation of the persistent currents in a 
normal loop 



B. Supercurrent in short SNS junctions 

In the case of short SNS junctions, the eigenvalue equa- 
tion Eq. (|44|) can be inverted. For the j = 1, we obtain 
e = Acos27tj; in the interval —1/2 < y < —1/4 and 
for j = —1, we have e = — Acos27ry in the interval 
< y < 1/4. These functions e>{y) are shown in Fig. 
IHt Now, we detail the calculation of the current for the 




d N « {„ 
ds » £„ 



FIG. 16: e J (y) for j = +1 (solid line) and j = —1 (dashed 
line) ds = oo, djv = 0. 



j = 1 case. From Eq. (33), we get 



„- =1 4A 1 7T777 f- 1/4 , 

l J m = cos h 27r / ay cos zirmy cos 27tt/ 

4> m y 2 y-i/2 

(69) 

Odd harmonics with m > 3 are zero. Even harmonics 
m = 2p are: 



^ = ^(-1)- 2P 



'2p 



4p 2 - 1 



(70) 



In this appendix, we show how the method developed 
in section III to calculate the harmonics of the current 



can be used to obtain the current of the purely normal 
ring. In this case, it is easier the equilibrium single parti- 
cle states spectrum rather than the excitation spectrum. 
For a free electron in a ring of length L, the electronic 
levels are: 



£n(^) 



2ir 2 h 



2*2 



elL 2 



(72) 



The derivation of Eq. (t 
rium spectrum. 



|) is still valid for the equilib- 



Irn(T = 0) 



4 1 

7rm 4> 



de 
dy 



(yp) cos 2-Kmyp 



dy-p^— cos 2irmy 



d 2 y 



(73) 



We have added an additional factor 2 for the spin de- 
generacy. The first term in Eq. @ gives the above 
result Eq. ( f43| ) but the second term give a non vanish- 
ing correction due to the parabolic dispersion relation. 
This correction is easy to evaluate since the curvature 
d 2 e/d 2 y = 4:TT 2 h 2 /m e iL 2 is a constant and we obtain: 



4 evp 1 v sinmfcpi 
61 </> = TTTT E i—sin2Trm<t> <f> 74 

7T L KwL * — ' 



m—1 



which is the l/fc^i-corrective term to the th order term 
obtained in Eq. (ph . 
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